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Resummed one— loop determination of the phase boundary of the SU{3)ji x SU{3)l linear 
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T. HerpajQ 

Department of Physics of Complex Systems, Eotvos University, H-llll Budapest, Hungary 

Zs. SzefQ 
Research Institute for Solid State Physics and Optics of the Hungarian Academy of Sciences, H-1525 Budapest, Hungary 

Complete one-loop parametrization of the linear sigma model is performed and the phase boundary 
between first order and crossover transition regions of the m,^ — mx-plane is determined using the 
optimized perturbation theory as a resummation tool of perturbative series. Away from the physical 
point the parameters of the model were determined by making use of chiral perturbation theory. 
Along the diagonal m^ = mx of the mass-plane we estimate m^ = 110 ± 20 MeV. The location of 
the tricritical point on the m^ = axis is estimated in the interval rn^ £ (1700, 1850) MeV. 
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INTRODUCTION 



In an attempt to understand the restoration of chiral and axial U{1) symmetries, chiral effective models are actively 
investigated (see e.g. [l|,y,y| for some recent works). Effective models indicate a very rich structure for the strongly 
interacting matter as function of quark masses and various chemical potentials [j,|^. The effective treatment represents 
a complementary approach to the lattice QCD field theory which, however based on first principles, has difficulties 
^ . mainly related to the computational power, in going towards the chiral Hmit to„ = m^ = m^ = 0. These effective models 
are constructed to share the same global symmetries as the massless QCD. It is expected that the lower m„, m^, Wg 
quark masses are (or alternatively m.,^ and ttik) the better they work. Universal arguments ja| predict a first order 
phase transition for the chiral limit. Lattice simulations with staggered quarks with a pion to rho mass ratio tuned to 
its physical value demonstrate a crossover type transition [3|- 

In QCD the critical line separating first order transitions from the crossover region in the m„ ^ — ms-plane is not 
precisely mapped, because of the difficulties of simulating dynamical fermions. There are several lattice studies with 
degenerate quarks to„ = md = tis , which show that the value of the pion mass on the boundary between the crossover 
and first order phase transitions drops substantially when finer lattices and improved actions are used, from the initial 
estimates of m'^ « 290 MeV p| or m,^ « 270 MeV J3| to m'^ = 67(18) MeV [lO|- In view of such low values one hopes 
that the boundary of the phase transition can be investigated reliably using effective chiral models. 

Although in principle it is simpler to solve an effective model than QCD, an exact solution cannot be given. Finding a 
good parametrization and an adequate method of approximation are the key issues when dealing with them. Attempts 
to parametrize physically the linear sigma model (LaM) date back to the early 70 's when in a series of papers Haymaker 
and collaborators have performed it at tree-level and started to calculate one-loop corrections at zero temperature (see 
[ll| and references therein). Recently other parametrizations were proposed in the literature [l2, UjI (see also [l3|l. 
It turned out that at tree-level it is not possible to fix the parametrization of the model using only the well-known 
pseudoscalar masses, information is also needed from the less known scalar sector. Moreover, the consequence of 
performing a tree- level parametrization is that one omits the effect of zero temperature vacuum fluctuations, which 
logarithmically depend on the renormahzation scale. At flnite temperature in the broken symmetry phase, the omitted 
terms have an additional implicit dependence on the temperature through the masses which depends on the order 
parameter. If the effective model is solved in an approximation which is not renormahzation scale invariant, then 
the renormahzation scale appears as any other parameter of the theory, and it has to be included in the process 
of parametrization in which some quantities calculated at quantum level are matched against their experimentally 
measured physical values. The effect of the renormahzation scale turns out to be both quantitatively and qualitatively 
important. It can have an effect on the pole structure of the scalar Green's function in the complex plane, as it 
happened in Ref. [I^. It influences the temperature dependence of the vacuum expectation value, and it can happen 
that above some temperature there is no solution to the equation of state (see e. g. [la|)- The renormahzation scale 
can even change the order of the phase transition. All this reflects the approximate nature of the solution. A good idea 
is to choose a range of the renormahzation scale where its variation affects the other parameters of the theory and the 
physical quantities less (e.g. trying to achieve approximate renormahzation scale independence). 

Due to the effects of the renormahzation scale mentioned above, it seems customary in the literature not to use the 
zero temperature quantum fluctuations of the fleld theory, and forgetting the renormahzation issue just solve the model 
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in a statistical mechanics inspired finite temperature quasi-particle approximation. Attempts to renormalize the model 
in the Hartree approximation of the CJT-formalism [13 were reported in [l^, but the result was not satisfactory. 
Recently much effort has been put in the renormalization of self-consistent resummation schemes of finite temperature 
QFT [13, I23, I21I I23I ■ In view of these results, solving the model in a properly renormalized approximation is nowadays 
a compelling requirement. In the present paper we want to go beyond the tree- level treatment of the model and its 
quasi-particle thermodynamics as it was treated in [lj|, and investigate the challenges of solving the renormalized 
version of the model by taking into account the logarithmic corrections. In particular, we want to investigate the extent 
they infiuence the location of the phase boundary in the pion-kaon mass-plane. 

In section ini we present the one-loop parametrization of the model in the ?7i,r — TO^-plane. It turns out to be rather 
hard to find a unique parametrization which works in the relevant region. The thermodynamics and the infiuence of 
the logarithmic terms are discussed at the physical point in section UTTl In section |^ we describe our results on the 
phase boundary, and we conclude in section Fvl 

II. PARAMETRIZATION OF THE MODEL AT ONE-LOOP LEVEL 

The Lagrangian of the ^[/^(S) x SUr{3) symmetric linear sigma model with explicit symmetry breaking terms is 
given by 

L{M) = ^Tr [d^^M^ni + iji'M^M) - h ('R {MHl)f - /2TV [M^ Mf - g (det(M) + det(Aft)) + e,,a., + ^yOy, (1) 

where the mixing sector is written in the non-strange (x)- strange (y) basis instead of the original 0-8 basis, by 
performing an orthogonal transformation on the fields as in [ij| (see AppendixEl- The complex 3x3 matrix M defined 
by the scalar (cr) and pseudoscalar (tt) fields can be written as 

1 ^ 1 

M = — = '^^{cfi + i7r,)A, + — =diag(<7a; + ii^^, a^ + itTx, \/2{ay + iny)), (2) 

^ ?'— 1 ^ 

where Xi : i = 1 . . .7 are the Gell-Mann matrices. Isospin breaking is not considered, therefore in the broken phase 
only the scalar fields a^ and ay have non-zero expectation values: x := {<Tx), y ■— {^y). After shifting the fields in the 
Lag rangian by their expectation values with a little bit of algebra one can perform the traces. Details can be found in 
[llirUfl|. Requiring that the sum of terms linear in the fiuctuations vanishes we obtain two equations of state. They 
are given explicitly in section ITTTI The coefficients of the quadratic terms are the tree- level masses (see Tab.^l, while 
the third and fourth order terms give the three- and four-point interaction vertices. 

In what follows, a set of non-linear one-loop equations will be given which determines at T = the 8 parameters of 
the Lagrangian: the couplings ^, /i, /?, q, the condensates x, y and the external fields e^, Cj,. Many ways of selecting 
these equations can be envisaged, see [ll| for alternatives. We have chosen to use as input the low lying pseudoscalar 
mass spectrum, namely the pion, kaon and eta meson masses and the decay constants of the pion and kaon, because 
they are the best known theoretically. 

In the broken phase a resummation is need/home/szepzs/fermion/renorm/kor3/proof/DS9047-l.psed, in order to 
avoid the appearance of negative mass squares in the finite temperature calculations of one-loop quantities. This can 
be done for instance using the Optimized Perturbation Theory (OPT) of Chiku and Hatsuda [1^1 . In the OPT the 
mass parameter — ^^ of the Lagrangian, which in the broken phase could be negative, is replaced with an effective 
(temperature-dependent) mass parameter m^ which is determined using the criterion of fastest apparent convergence 
(FAC). The mass term of the Lagrangian reads: 

Lmass = -m^TrM^AI - -{^i^ + rr?)T:v mHi = -m^TxMHl - -Am^Tr AfU'/, (3) 

where the finite counterterm Am^ is taken into account first at one-loop level. 

This resummation method replaces — ^^ by the effective mass square rr? in the tree-level masses (see Tab.H)), and 
preserves all the perturbative relations upon which Goldstone's theorem relies [la|- The renormalization is achieved 



ml =m^ + 2(2fj+f2)x^+4fry^+2gy 


ml^ =m2 + 2(2/1 +3/2)x2 + 4/iy2-2py 


ml = m^ + 2(2/1 + h){x^ + y') + 2./-22/' - ^/2x{2hy - g) 


ml = m^ + 2(2/1 + f2){x^ + y^) + 2/2J/' + ^/2x{2f2y - g) 


mt. =m' + 2(2/1 + f2)x^ + 4/iy2 _ 2gy 


"^Lx = m^ -f 6(2/1 + f2)x^ + 4/1 y2 + 2gy 


ml,,=m^+4f^x^ + 4{,h + f2)y^ 


ml^y = m^ + 4/1 a;2 + 12(/i + f2)y^ 


^Ly = -25a; 


"T-Lv = ^h^y + 25a; 



TABLE I: The squared masses of the (pseudo)scalar nonet appear in the (first) second column. The last three rows represent 
the mixing sectors. They can be written in the conventional basis using the formulas of Appendix 1X1 



both in the symmetric and the broken phase by the following counterterms 



13A^ + 12/1/2 + 3/1 , A2 



3/1/2 + 3/1, A2 
22 



ej. «JJ1J2 -r 0^2 1 ^'^ ,r,\ 

5f2 = 3 In-^, (7) 



where A is the 3d regularization cutoff, and / is the renormalization scale. Note that only the mass counterterm differs 
from its standard expression [ll|. In the present form this counterterm is temperature-dependent through the effective 
mass, but this temperature dependence is canceled by higher-loop terms |lal23|- ^^ what follows, all quantities and 
equations are renormalized without any change in the notations. 

The above mentioned FAC criterion, which determines the effective mass square m^, is realized in the present case 
by the requirement that the pole and the residue of the one-loop pion propagator 

stay equal to their tree-level values. Here we anticipated that we also need a finite wave function renormalization in 

_i 

order to make the residuum equal to 1, and rescaled the pion fields as tt — *■ Ztt ^ tt. 

According to this FAC criterion the inverse of the finite wave function renormalization constant is 



^-1 .^ ^ _ as^(p2,TO,,o 



(9) 



The one-loop pion pole mass 

Af2 = -^2 ^ (4^^ ^ 2/2)x-2 + 4/1J/2 + 2gy + Re {^^(p^ = Mlm.J)} (10) 

I 
has to be equal to its tree-level value {Mt^ = rriT^). Therefore, using the expression of the tree-level pion mass of Tab. HI 

the following "gap" equation can be obtained for the effective mass: 



m 



-m' + Re {^Ap^ = ml,m,{m^),l)} , (11) 



where the m^-dependence of the self-energy (through the tree-level masses) is explicitly shown. The different contribu- 
tions to the self-energy are depicted in Fig. Q] 

The effective mass can be replaced by the pion mass by expressing it from its tree-level formula. Then ifTTll can be 
interpreted as a zero temperature gap-equation for the pion mass: 

ml = -^i^ + (4/1 + 2/2)a;2 + if^y^ + 2gy + Re {^^(p^ - ml,m,imA, I)} , (12) 

where the tree masses of all mesons are expressed through the pion mass. A similar gap-equation will be used in the 
thermodynamical calculations for the temperature dependence of the pion mass. At T = 0, the task is "reversed": the 
pion mass is known and (I12|l belongs to the set of equations, which determines the parameters. We have chosen to 
express the effective mass m^ from the tree-level mass formula of the pion because the pion has the smallest mass, 
and positive solutions of ltT2)l ensure the positiveness of all the other masses. We use the kaon and eta masses, the 
relations of the Partially Conserved Axial- Vector Current (PCAC) for the pion and kaon at one-loop order, and the 
two equations of state to fix the remaining parameters. 
The one-loop kaon propagator is the following: 

,■7-1 
Z]^^ and the one-loop pole mass of the kaon Mk can be calculated similarly as in the case of the pion: 



Z^^ 



-1 



dT.K{p ,mi,l) 



and 



(14) 



o^^M? 



Ml = -^? + 2(2/1 + h){x' + y2) + 2/2y2 - V2x{2hy - g) + Re {S^l/ = Ml,m,, I)] . (15) 
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FIG. 1: The physical content of the one-loop pseudoscalar self-energies. 



The description of the t] and rj' mesons is sUghtly more comphcated because of the mixing in the x — y sector (0 — 8 
in the conventional basis). The propagator is a 2x2 matrix, and pole masses are defined as the real part of the solutions 
of the following equations 



Det 



V - m'^,, - E,,^^ {p^,m„ I) -m\^^ - S,,^^ (p^ m„ I) 



-m, - 



p2^A/2,M^, 



^ri^yiP ,mi,l) p^ -m^^^^ - ^,j^^{p^,mi,l) 
This yields two equations for the mass eigenvalues M^, Af^-: 

M2 = iRe [ml^ + S^^^ (p^ = M^ m„ /) + m^^ + E,^„ {p^ = M^ m„ /) 



(16) 



A/ 2 



-0)^},(17) 






^Re {m^^^ + E^^^ [p" = M^^,, m„ + m^^^ + E„^„ (p^ = Af^ ^ m„ /) 



+ ^{m^^^+^^.AP^ = Af 2, ,„,,/) _™2^^ _ S,^jp2 ^ M2 ,m„0)2 +4(m2^^ +E,^„(p2 ^ A/^ , ,„,, ;))'}(-18) 

The definitions lll()|l , (jl5|l , Ijl7|l and (I18|l of the pole masses give the correct one- loop masses only when the self-energy 
is not complex. We note that if the tree-level masses are close to their experimental values, then by looking at the 
self-energy contributions of pion and kaon in Fig. ^ one can recognize that they have no imaginary part for p^ = m^ 
and p^ = m\ , respectively. This is also true in the case of the ry self-energy. It turns out that the rj' self-energy has 
an imaginary part at the pole-mass determined as the zero of the real part of the self-energy except the narrow range 
of 1820 MeV < I < 1880 MeV. For this reason we decided not to include the one-loop equation for M,f into the set 
of equations used for the parametrization. We make up for the missing equation by extending FAC criterion to the 
one-loop kaon mass too. This condition reads: 



a4 = m|- = ml~ 2gy + 4/2^^ - V2x{2f2y - g). 
Two more equations are provided by the one-loop PCAC relations which according to [T^l reads 

-^-iD-^(p^Q) 



as 



M2 



Ik 



-^-iDy/{p^0)x + V2y 



'K 



'Ji 

Ml 



(19) 

(20) 
(21) 



As shown in Appendix El these equations can be rewritten in an explicitly renormalization scale- independent form. 
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FIG. 2: The renormalization scale dependence of various quantities at the physical point: the parameters (a), the non-strange 
(x) and strange (y) vacuum expectation values and the finite wave function renormalization constants Z^v and Zk (b), the 
pseudoscalar masses (c), and the scalar masses (d). 

Finally, the last two parameters, the two symmetry breaking external fields e^ and e^ are determined by the one-loop 
equations of state, with the help of zero temperature chiral Ward identities (Appendix [B)|: 



e, = Z-i(-zi?-i(p = 0)); 



'K 



{-iD-^\p = 0)) (^-^+y^- Z-' {-iD-\p - 0)) 



X 



(22) 
(23) 



One can notice that, since OPT preserves Ward identities at tree and at one-loop level as well, the above parametriza- 
tion, in which the tree-level masses of pion and kaon equal the one-loop masses, ensures at zero temperature the validity 
of Goldstone's theorem for both pion and kaon. 



A. Parametrization at the physical point 

The parameters were determined as follows. From (IIOII and IjlOp one can express /i^ and g, respectively. Next, from 
the system of 4 non-linear equation (tlfil . l(T7j) . lf2(l . (l2T|) one can numerically determine /i, /2, x, and y as functions 
of the renormalization scale I. Going back to ifTIijI and ifToll one can compute y?' and g, respectively. Substituting these 
parameters into Ij22|l and \'2?>\ one can determine €x and Cy. Unlike the tree- level parametrization case [lj|, now all 
scalar masses are predicted. 

The numerical solution for different renormalization scales I is presented in Fig. [21 (a) for the physical point, where 
TO^ = 138 MeV, rriK = 495.6 MeV, m^ = 547.8 MeV, /^ = 93 MeV, and Jk = 113 MeV. In figure Fig.|2(b) one can 
see the renormalization scale {I) dependence of the non-strange {x) and strange (y) vacuum expectation values and 
of the finite wave function renormalization constants Z-^ and Zk- They have a plateau for I < 1400 MeV, and the 
tree-level m^' (see Fig. [21(c)) is the closest to its physical value in the region I g (1000, 1400) MeV, where the variation 
of the tree- level scalar masses (see Fig. [3 (d)) is the mildest too. We have decided to use in our thermodynamical 
investigation this range of the renormalization scale in which the tree-level masses entering into the propagators of 
Fig- Hare reasonably close to their experimentally measured values. In Fig. [21 (c, d) we present an estimation of the 
predicted one-loop pole-masses based on the real parts of the corresponding self-energies. This is a good approximation 
in the case of rj' , uq and /o since the zeros of the inverse propagators correlate well with the location of the well-defined 
peak in the corresponding spectral functions. The one-loop mass Mfg, which is not shown in the figure, has a rather 
large value in the present range of the renormalization scale (decreasing from Mf„ = 2000 MeV for / = 1000 MeV to 
Af/o = 1400 MeV for I ~ 1400 MeV). The shapes of the spectral functions of k and a (see Fig. [21 for pa) are more 
complicated, they have a threshold dominated peak with large width, and are very sensitive to the renormalization 
scale. In this case, it would be more appropriate to define the mass and width of a decaying particle as the real and 
imaginary part of a complex pole. In the 0{N) model in the large N approximation J2J|, this continuation into the 
second Riemann sheet was performed in the sigma channel and the poles of the propagators were determined. In this 
model the continuation of the propagators into the complex plane would be more difficult due to the appearance of 
many decay thresholds and is beyond the scope of the present investigation. 
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FIG. 3: The kaon mass dependence of m,, as given by ChPT and the Veneziano formula for m^n- = 10 MeV. 

B. Parametrization in the m-^ — ttik plane 

Since we are interested in the phase boundary on the ititt ~ m/f-plane, we have to take into account the variation 
of the parameters with TOtt and rriK- A method for the parametrization away from the ph ysic al point was proposed 
in [iJI, which rehes on the formulas provided by the Chiral Perturbation Theory (ChPT) [2^. Because our present 
parametrization does not use the rj' meson, we make use of the SU{3) ChPT describing the chiral dynamics of the 
pseudoscalar octet. In the large Nc limit, the formulas for the pion, kaon mass dependence of the decay constants and 
of the 77 mass up to 0{1/P) read as |2fl |: 



fK 



m^ 






4to|. 



f(2L. 



L. 



'K 



-Kf 



P 



(24) 
(25) 

(26) 



where L5 and Lg are low energy constants and / is the decay constant in the chiral limit. All the parameters of the 
large Nc limit of the SU{2>) ChPT can be determined at the physical point from the equations above. Their values, 
X5 = 2.0152 • 10-3, Lg = 8.472 • 10"^ and / = 91.32 MeV, are fixed for all values of m^ and m^. 

For low values of m^ and niK the sensitivity to the renormalization scale of the LaM is bigger than the uncertainties 
coming from the omission of the chiral logarithms. As an effect of these chiral logarithms, for large values of rriK, the 
formula of the SU{3) ChPT yields a decreasing value for to,, for increasing ttik, see Fig. |3| If we use lf26jl . the same 
behavior occurs at a larger value, which is around uik ~ 1300 MeV. This is non-physical as both the kaon and the eta 
particles have to decouple in order to arrive at the 0(4) model for rriK -^ 00. This shows the failure of the ChPT at 
high values of the kaon mass. 

In view of the bad behavior of to,, determined using large A^c ChPT, we have used, as an alternative, the following 
mass-formula by Veneziano [23 : 



2 2 

TO^ =TO^ 



2^^vO 




2ml 



g^'^^O 



9^Ko- 



(27) 



Ato^q is the non-perturbative gluonic mass contribution in the singlet channel of the mixing i] — rj' sector, related to 
the axial C/(l) dynamics. Using the values of the masses at the physical point in lf27jl one can fix the value of the extra 
mass contribution: Am^g = 2.3 GeV^. One can see in Fig. El that this parametrization gives for to,, values which are 
almost identical to the values coming from the formula of ChPT in the large Nc limit, up to values of uik for which 
ChPT breaks down. 

We note, that in the original paper Ato^q was determined using the trace of the 2x2 matrix of the mixing rj — -q' 
sector. We indulged in modifying the procedure in order to make contact with the ChPT in the large Nc limit, as to,, 
obtained form 112 7|l with the original parametrization is always smaller than the value given by the large Nc ChPT. 

The continuation onto the to^ — toj^ -plane of fyr, /k, and to,,, based on the formulas of this subsection, allow us to 
determine the parameters as described in III Al in a wide region of the mass-plane, except for high values of niK, near 
the m^ = axis. 



III. COMPLETE 1 LOOP THERMODYNAMICS OF THE LcrM 



With the intention of determining the order of the phase transition in the pion-kaon mass-plane we have to monitor 
the order parameters as functions of temperature. They are obtained from a set of three equations: two equations of 
state for x and y and a gap-equation for the pion mass. The temperature dependence of the order parameters at finite 
T can be obtained from the equations of state: 



e^ + m^x + 2gxy + Afixy'^ + (4/i + 2f2)x^ + ^ ,^^1(1, m,{T), T) + llm^x = , 

i 

-ey + m^y + gx' + Ahx^y + A{h + f2)y^ + Y,J^tpil,MT),T) + Am^y = 0, 



(28) 
(29) 



where in the tadpole integral I{l,mi(T),T) we have explicitly displayed the implicit temperature dependence of the 
tree- level masses, expressed with the pion mass determined by the gap-equation 



= -^i^ + (4/i + 2h)x^ + 4/iZj2 + 2gy + E,(p2 = ml,m,{m^),l) + S^(w = 0, m,{m^)). 



(30) 



The sum goes over all mass eigenstate meson fields with isospin multiplicity factor Ji: J,r.ao =3, Jk,k = 4, and 
Jri,7j',ajo — 1- The coefficients tf and t^ appearing in ll28)l and lf29jl are listed in Appendix C of [lj|- The standard 
one- loop integrals appearing in the formulas above can be found for instance in [ia|. In the present study we use the 
one- loop bubble integrals appearing in (j3()ll for |p| = 0, corresponding to particles at rest. As in zero temperature 
case, at finite temperature OPT guarantees through the gap-equation the validity of Goldstone's theorem for the pion. 
Following Ref. [io|, S^, the finite-temperature part of the self-energy, is taken not on the mass-shell, but instead at 
w = 0. This is done because above a given temperature S]^(r)(a; — m.,^{T), 0) becomes complex, and the real solution 
of the gap equation ceases to exist. At the physical mass-point this temperature is typically below the value of the 
pseudocritical temperature, invalidating the study of the phase transition. The imaginary part is produced by one-loop 
bubble integrals where two unequal masses mi and 771,2 appear, when the relation uj = m-r; < \mi — m2\ is satisfied. 
It has the consequence that we can not regard the pion gap-equation as the equation determining the "true" one-loop 
mass of the pion. Due to the imaginary part of the self-energy the most adequate way of determining the pion mass 
would be to look for a complex pole of the pion propagator. Still, in the present work we content ourselves with the 
study of the spectral function of the pion, which, as we will see, in certain temperature ranges also provides information 
on the pion mass. 



A. Influence of the logarithmic terms 



In order to estimate the infiuence of the logarithmic terms on the solution of the equations of state we have performed 
the following check. We have taken the parameters determined with the one-loop parametrization presented in Sec.lTTl 
and modified the value of //^ such as to incorporate the zero temperature logarithmic terms. Then, at finite temperature, 
we have solved the gap equation and the two equations of state without taking into account the logarithmic terms, which 
implicitly depend on the temperature through the masses. The difference between these expectation values and the ones 
calculated with the logarithmic terms (and with the original, unmodified value of ^^) is only due to the temperature 
dependence of the logarithmic terms. This can be seen in Fig. ^ without the logarithmic terms the pseudocritical 
temperature is lower by about 20%. The variation of strange condensate with the temperature is significantly different 
in the two cases. We notice that at high temperature the solution obtained with logarithmic term included ceases 
to exist, as was also observed in [l^. This is a non-physical phenomenon, it happens before the restoration of chiral 
symmetry completes. We consider the solution reliable up to temperature values which are below the turning point in 
X and y. 
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FIG. 4: Comparison of the temperature dependence of strange (y) and non-strange (x) condensates with and without the 
inclusion of the logarithms for I = 1200 MeV. 
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FIG. 5: The temperature dependence of non-strange (x) and strange (y) order parameters for different values of the renormal- 
ization scale (l.h.s) and the T-dependence of the tree-level masses for I = 1200 MeV (r.h.s). 

B. The solution at the physical point 

At the physical point and in the investigated range of the renormalization scale, the behavior of the non-strange 
order parameter shows a smooth restoration of the SU{2) chiral symmetry. The pseudocritical temperature moderately 
depends on the renormalization scale. The strange order parameter varies less and it is approximately scale-independent 
until the solution is reliable, see the left panel of Fig.0 The right panel shows that the tree-level SU(2) mass partners 
tend towards degeneracy as the temperature increases. Unfortunately the solution falls dead before the restoration of 
the complete SU{3) symmetry, due to the effect of the logarithmic terms f see IIII A)l . 

Fig. El illustrates the temperature dependence of the spectral functions in the pseudoscalar and scalar channels and 
the behavior of the zeros of the real part of the inverse pion and sigma propagators. At low temperature the spectral 
function in the pion channel develops a peak whose location is close to the physical mass value of the pion. Close to 
the pseudocritical temperature there are significant changes in the peak structure. In this temperature range, based 
on the spectral function, one can not determine the true pion pole-mass. In this range, the tree-level pion mass, which 
is the solution of the gap-equation lfT2|l interpolates between various zeros of the inverse propagator, see the left panel 
of Fig. El At high temperature the spectral function has a well defined peak again, whose location correlates with the 
zero of the inverse pion propagator. 

At low temperature the spectral function in the sigma channel has a large width and the zeroes of the real part 
of the inverse propagator are very sensitive to the renormalization scale, therefore the sigma pole mass can not be 
estimated. We can define the pole mass only at high temperature, where the spectral function develops a well-defined 
peak, correlated with the zero of the inverse sigma propagator. Increasing the temperature its location approaches the 
location of the peak in the pion channel and eventually they become degenerate. This shows that the one-loop masses 
of the pion and sigma also reflect the restoration of the SU{2) chiral symmetry at high temperature. 
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FIG. 6: The spectral function of the pion (l.h.s.) and sigma (r.h.s.) for various values of the temperature and I — 1200 ]V[eV. 
The zeros of the real part of the inverse pion/sigma propagators and the corresponding tree-level masses are also depicted, as 
lines in the T — u> plane. 



IV. THE PHASE BOUNDARY ON THE m^ - m^ PLANE 

The main result of our work is the determination of the boundary between the region where a crossover transition 
occurs, with a smooth variation of the order parameters, and a first order phase transition region, which is signaled by 
the multivaluedness of the order parameters. 

It is remarkable, that at large values of the kaon mass, the boundary proves to be independent of the renormalization 
scale ^. For large values of uik the existence of a scaling region belonging to a tricritical point (TCP), with mean-field 
exponent can be confirmed. This is a new feature of the complete QFT treatment, it was not observed using tree-level 
parametrization! From mean-field studies (see e.g. J23|) it is known that near the TCP the boundary of the edge of the 
first order region deviates from the to„ ^ = axis of the quark mass-plane according to m.u^d ~ ("t-^*^^ — rris)^'^. Using 
the tree-level formulas of ChPT, for m^ and m^, namely rn^ = 2mi3o, rrij^ = (m + ms)Bo with m = ^{rriu + md), 

one can easily translate this into a relation between the critical values of rriK and 771,^: m^ = m'j^\^ + ^ — am J . 
Due to the failure of our parametrization at very high values of itik , and near the uIt^ = axis we use this formula to 
extrapolate the upper edge of the phase boundary to the Wif-axis. Using large Nc ChPT theory we obtain mjf'^ = 1718 
MeV, while using the Veneziano formula lf27ll we get mjf^ = 1838 MeV. Both fits are of very good quality suggesting 
the radius of the scaling region to be Ato^ « 40 MeV. In terms of the strange quark mass our estimate corresponds to 
jy{rcp = 13 _ 15 X 771s . This result can be compared to the recent lattice result of J23|, where it was estimated that 
jy^tnc j^ 2j^^ ^ jj^ pg^gg q£ J2^ ^YiQ piece of the phase boundary used in the extrapolation to rn^'^^ was much closer to 
the physical point and the location of the TCP was estimated using points with m^ < mPj^y^ . The lattice estimate for 
the location of TCP would improve if points could be simulated closer to the scaling region. 

In the present treatment the critical pseudoscalar mass for degenerate quarks was obtained in the range m!^ e (90, 130) 
MeV. The spread is the result of the renormalization scale dependence. This range of m!^ is significantly above the 
one previously obtained in [ij| using a tree-level parametrization, which eventually led to m% G (20, 60) MeV. It is 
interesting to note, that the boundary lies to the right from the fitted curves which describe the scaHng behavior near 
TCP. 
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FIG. 7: Phase boundary of LcrM obtained using a one-loop parametrization at T = and the formulas of large Nc ChPT for 
continuation on the m^ — TTiK-plane. For large values of ttik we have also used the Veneziano formula for tti,, cf. IIIBI 



^ It can be seen after some straightforward calculation, that as the kaon mass increases, the relative weight of the terms containing log(Z) 
are becoming negligible in the equations used for parametrization as well as in those used for thermodynamical calculations. 
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V. CONCLUSION 

In this paper we have studied the phase boundary of the Hnear sigma model in the m^ — ttik plane. We have 
managed to find in a resummed perturbation theory a set of renormalized equations based on which complete one-loop 
parametrization of the model is possible not only at the physical point, but in a large region of the mass-plane too. 
We allowed for the variation of the parameters with TOtt and ttik when moving away from the physical point, and used 
formulas of the ChPT for the continuation of certain physical quantities into the mass-plane. 

As a result of the one-loop solution of the model we were able to estimate the location of the tricritical point on the 
m^ = axis of the mass plane at a rather high value of the kaon mass. We have shown evidence of the existence of a 
scaling region around it. These were not observed previously in the quasi-particle approximations to the model. For 
the same values of rriK the phase boundary line lies at higher values of pion masses than in the case of the tree-level 
parametrization of the model, and without the inclusion of the vacuum fluctuations. For degenerate pseudoscalar 
masses the phase boundary was obtained in the range to^ e (90, 130) MeV. 

We have discussed also the limitation of the optimized perturbation theory, related to the fact that one looks for real 
solutions of the gap-equation for the resummed mass in temperature ranges where the self-energies develop imaginary 
parts too. It would be interesting to study the phase boundary using a self-consistent approximation based on the 
exact propagator, Hke the 2PI approximation. 
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APPENDIX A: CONNECTION BETWEEN (0-8) AND (X-Y) BASIS 

From the x — y basis used in section ^ one can obtain, with the help of an orthogonal transformation, the mixed 
scalar-, pseudoscalar and the external fields in the conventional (0 - 8) basis as: 

= 0(^" ) , whereO:=^ (^ U. (Al) 

'y/ v"»/ \"y/ V**/ \^y/ "^"^ \ ^ ^v2y 

The transformation of the mass matrices, and the mass eigenvalues are 

^2 





, m'i rni / I mi m^ ' ' ^ '^ 



+ ml ± ^ (mi - ml Y + Amf^ 

lc 'lyy y ^ 'Ixx 'lyy 'Ixy 



r;oo ^08 ^ (9 v.. "l^y ] O , and m^, ^ ^- . A2 

I ml / [ml ml ' ' " ■'' '^ ^ ' 

I?08 »?88/ \ Vxy riyy ^ 



Similar expressions also hold for the scalar fields. 



APPENDIX B: WARD IDENTITIES 

From l|Hl and ((T3|l . the inverse pion and kaon propagators at zero external momenta are 

7 — In — 1/"^ o,\ ,.^2 I v:^ / n\ I A. ^2 



^iZ-'D-'{p = Q) = mi + ^Ap = ^) + ^m' (Bl) 

^iZ]^^D-^\p = Q) = m|. + EA-(p = 0) + Am2, (B2) 

where the finite counterterms of OPT are now expHcitly indicated. Comparing the expressions of tree-level masses 
(Table nj with the tree-level part of the equations of state Ij28|l . (I29|l . one can obtain the corresponding tree-level Ward 
identities e,jr = mja; and e^ = TO^(-\/2/2.T-|-y) — m^-\/2/2a;. The diagrams for the self-energies are shown in Fig.^ They 
include both tadpole and bubble diagrams, and the latter ones can be decomposed for non-equal propagator masses 
and zero external momenta into the difference of two tadpoles: 



.Q. 



(B3) 



Therefore, the pseudoscalar self energies at p = can be represented as a linear combination of tadpole contributions 
whose weights are complicated expressions of the tree-level masses and the corresponding four- and three-point couplings 
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(see [U). In the case of pion and kaon these weights simphfy to 

E.b = 0) = T.J^tfI{m.,l)+.Am^ 

X 

^ , n^ E J^(M + ^t\)I{m.,^) + (x + 72y)Am^ 

^^-(^ - 0) ^771^ ' (B5) 

where I {mi, I) is the T = tadpole integral and the sum goes over all mass eigenstate meson fields with isospin 
multipHcity J^. Substituting these and the corresponding tree- level Ward identities into the expressions l|Hl, CSjl, 
one can obtain the equations of state lf28jl and lf29jl . which determine the external fields at zero temperature in the 
parametrization process. The relations IJ22II and ^2'^\ are also valid at finite temperature and ensure the fulfillment of 
Goldstone's theorem at one-loop order. 

APPENDIX C: PCAC RELATIONS 

The one-loop order PCAC relations for the pion and kaon are given in [T^l by 



UM^ = ^Z^e, , fKMf, =. ^-^ + ^Z^e, . (CI) 

With the help of Ij22|l . Ij2^^|l the external fields can be eliminated and one can obtain the expressions II2()|I . (I21|l . which 
appear renormalization scale-dependent. However, one can rearrange the pion self-energies appearing in (j2()ll as follows: 

_ ml + ^^ip = OJ) ^^^ _ ml + S^(p2 = Mil) + (S^p ^ 0, /) - S.(p^ = M^ I)) ^- 1 

^(l-^-^4^)^.- (C2) 

In the last step above we used the fact, that mj + St(p^ = ^^l is just the definition of the pole mass A'/^. The 
p-dependent part of the self energy, E, like the wave function renormalization constant, does not depend on the 
renormalization scale and in consequence ljC2ll is actually scale- independent. The PCAC equation for the kaon can be 
analyzed in a similar way. 
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